Statistics at the tip of a branching random walk and the delay of traveling waves 
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We study the limiting distribution of particles at the frontier of a branching random walk. The 
positions of these particles can be viewed as the lowest energies of a directed polymer in a random 
medium in the mean-field case. We show that the average distances between these leading particles 
can be computed as the delay of a traveling wave evolving according to the Fisher-KPP front 
equation. These average distances exhibit universal behaviors, different from those of the probability 
cascades studied recently in the context of mean field spin-glasses. 
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The interest for branching random walks has a long 
history in Mathematics [H [3 E] j Physics and Biology. In 
Biology they are commonly used to model the genealogies 
of evolving populations, the spread of an advantageous 
gene or of an infection, the combined effects of selection 
and mutations [U |5J [B]. In Physics they also appear 
in many contexts such as reaction-diffusion models [7, 
[8], particle physics [Qj [10], or the theory of disordered 
systems [111112]. 

In one dimension, the right frontier of a branching ran- 
dom walk is the region located near its rightmost parti- 
cle. An interesting question is what does the branching 
random walk look like when seen from this frontier. For 
example one can try to determine the position of the 
second, the third, ... or the n*^ rightmost particle in the 
frame of the first rightmost particle. The statistical prop- 
erties of these positions depend on time and have a well 
defined long time limit [5] which we study in this letter 
using traveling wave equations of the Fisher-KPP type 

miaiii 



dh 
dt 



(1) 



The fluctuating distances between these rightmost par- 
ticles allows one to understand why directed polymers in 
a random medium |llj have non-selfaveraging properties 
similar to mean field spin-glasses [13]. Their study is 
also motivated by the growing interest for the statistics 
of extreme events [12lll6l[lI||l8lll9l[2Ql|2ll[22] which 
dominate a number of physical processes [231 121] • The 
last two decades have seen the emergence of universal sta- 
tistical properties of the probability cascades describing 
the energies of the low lying states of several spin-glass 
models [ 22] [25l [26l [SZllMl [MllSO] . Somewhat surprins- 
ingly, as shown below, the distribution of the distances 
between the extreme positions of particles in a branch- 
ing random walk (which are nothing but the energies of 
the low lying states in the mean field version of directed 
polymer problem [llH is different from the predictions of 
the probability cascades [22, 25j |26l EH [28 , 29 , 30 . 

To start with a simple case, we consider a continu- 
ous time branching random walk in one dimension: at 



time t — Q there is a single particle at the origin a; = 0. 
This particle diffuses (for convenience we normalize the 
variance of its displacement during time t to be 2t) and 
branches at rate 1. Whenever a branching event occurs, 
the offspring become themselves branching random walks 
which diffuse and branch with the same rates. The num- 
ber of particles grows exponentially with time and they 
occupy a region which grows linearly with time [T| [2] . 

It has been known for a long time [H |2l [12] that the 
distribution of the rightmost particle of a branching ran- 
dom walk can be determined by solving a traveling wave 
equation: the probability Qo{x,t) that, at time t, there 
is no particle at the right of of a;, satisfies 



dQo _ d^Qo 2 n 



(2) 



[The derivation of ^ is standard: one decomposes the 
time interval {0,t + dt) into two intervals (0,di) and 
{dt,t + dt), and write that Qo{x,t + dt) = Qo{x,t)^dt + 
{Qo{x — ri,t))n{l~ dt) where the first term represents the 
contribution of a branching event and 77 in the second 
term the displacement due to diffusion during the first 
time interval (0, dt). With our normalization {rf) — 2dt] 
Up to the change h = 1— Qo, equation Q is the Fisher- 
KPP equation ([ij. Since at time i = there is a single 
particle at the origin, the initial condition is simply 

Qo{x, 0) = 1 for a; > 0, Qq{x, 0) = for x < 0. (3) 

If Qnix,t) is the probability that there are exactly n 
particles on the right of x, one can see, as for Qq, that 
the generating function ipxi^x^t), defined as 



i>x{x,t) = ^A"Q„(x,t), 



(4) 



n>0 



evolves according to the same equation ([2|, the only dif- 
ference being that the initial condition is replaced by 

i^x{x, 0) = 1 for a; > 0, %l;x{x, 0) = A for a; < 0. (5) 

We are now going to see that the knowledge of ^\{x,t) 
allows one to obtain the average distances between the 
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rightmost particles of the system. If p„(a;,t) is the prob- 
ability of finding the n^^ rightmost particle at position x, 
it is easy to see that 

^^^pi{x,t) and = p„+i{x,t) - pn{x,t). (6) 

The average position {Xn{t)) of the n^^ rightmost parti- 
cle and the average distance {dn^n+i{t)) between the n**^ 
and {n+ 1)"^ rightmost particles can then be defined by 



{Xn{t)) = / XPn{x,t) dx, 



(dn,„+l(t)) = (^„(t))-(^„+lW). 



(7) 
(8) 



[One should notice that the normalization of pn{x,t) is 
not 1 but J Pn{x^ t)dx = (1— e~*)"~^ due to the events for 
which the total number of particles at time t is still less 
than n. One could prefer to use different definitions of the 
positions or of the distances, for example by conditioning 
on the fact that there are at least n -I- 1 particles in the 
system, but any such definition would coincide with ( 7|8 ) 
up to contributions which decay exponentially with time 
and disappear in the long time limit that we study below.] 
With the definition dSl we obtain from ( |4p|7t that 



E 



A"(d„,„+i(t)) = / X 



dQo _ dtpx 
dx dx 



dx, (9) 



which relates the distances {dn,n+i{t)) between the right- 
most particles to the solution ipxit) of the partial differ- 
ential equation ^ with the initial condition ([s]). 

We have integrated numerically the equations satis- 
fied by 'ip\{x,t) and its derivatives with respect to A to 
measure the distances (rf„,„+i(i)) between the n*'' and 
{n + 1)*^ rightmost particles. In our numerical integra- 
tion, we had to discretize space and time; we checked 
that our results shown in figure [l] were stable when we 
decreased our integration steps. 
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FIG. 1: The average distances between the first rightmost 
particles (di,2(t)), (rf2,3(i)) and (d3,4(i)) of a branching ran- 
dom walk versus for t up to 3000. 

One can first remark that, in contrast to standard 
Monte-Carlo simulations, where all the branching events 
would be simulated and for which the maximum reach- 
able time would be t '--^ 20 (with a number of particles e* 



of order 10^), the integration of ([2| or of its derivatives 
allows one to achieve much larger times. One can also 
notice in figure [l] that the distances converge like 1/i to 
well defined values. We will see that this Xjt convergence 
is consistent with our analytic expression ( |24[ ) below. 

We did not find an analytic theory to predict the lim- 
iting values that we measured as in figure [l] 



(di,2) 0.496, (d2,3) - 0.303, (^3,4) 0.219, 
(d4,5> - 0.172, (d5,6) - 0.142, (dej) - 0.121. 



(10) 



As shown below (17), we can however predict their large 
n behavior. 

Before doing so, it is interesting to compare our results 



( 10 1 to the expected values of the gaps between the low 
lying energies of spin-glass models such as the REM and 
the GREM [3TJ [32 • In these models one can show that 
these energies are given by probability cascades [22l [25j 
B71 [53] and that the energy gaps at the leading edge 
are the same as those of a Poisson process on the line 
with an exponential density. For such a process, with 
density e~"^, the probability distribution of the positions 
is Pnix) = exp[—nax — e~""^ /a] from which one gets 
through 



+ 1/GREM 



(11) 



Clearly there is no choice of a for which our numerical 



results (10 1 are compatible with (11). 

It is well known 1,2, J4^ that the solution Qq{x, t) of the 
Fisher-KPP equation (|2| with the step initial condition 
([3]) becomes, for large t, a traveling wave of the form 



Qo{x,t)^F[x-{Xi{t))], 



(12) 



where the shape F{z) of the front [F{z) \ as z ^ 00 
and F{z) ^ Q as z ^ —00] is time-independent and its 
position, which can be defined as the average position 
{Xi{t)) of the rightmost particle, has the following long 
time behavior p iTill55] 



(Xi(i)) =2t~ -lnt + 0(l). 



(13) 



ip\{x,t) is also the solution of the Fisher-KPP equation 
([2]), but with the initial condition ([5|. As ij\{x, 0) decays 
fast enough [M], one expects the same large t behavior 
as ( |12|13 1, up to a A-dependent delay /(A) due to the 
change of initial condition: 



Mx,t)^F[x^{Xi{t)) + f{\)]. 



(14) 



From ( |9|12|14 l we see that the translation /(A) is nothing 
but the generating function of the average distances 



/(A)- lim VA"(d, 



'n,n+l 



it))- 



(15) 



n>l 
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We were not able to find an analytic expression of the 
delay function /(A) for arbitrary A. For A close to 1, 
however, we are going to show that 

/(A) = TA-lnTA + 0(l) withTA = -ln(l-A). (16) 



This implies (151 that the distances have the following 



large n asymptotics 

t—oo 



(17) 



Compared with (jllj), we see that there is a correction, 
which we believe to be universal as discussed below. 
(Note that the same asymptotic distances would be ob- 
tained for uncorrelated particles distributed according to 
a Poisson point process with a density — a;e~^ for nega- 
tive X.) 



For A close to 1, the time t\ in ( 16 1 is the characteristic 



time it takes '4>\{—oo,t) to reach a value close to 0. The 



most nai've idea to derive (16 1 would be to say that it 



takes this time t\ for ipx{x,t) to look like the step func- 
tion Qq(x,0), and then to start moving like Qo{x, t). As 
the asymptotic velocity is 2 (see (13l) this would lead to 
a delay /(A) ~ 2t\ which is wrong (see (16l) by a fac- 
tor 2. The problem with this idea is that, while il}\{x,t) 
evolves to approach on the negative tail builds 

up on the positive x axis which has a strong influence on 
the dynamics later on. 

To derive ( [l6| , we need to understand the shape of 
^\{x,t) for t > T\. Let Yt be the position where 
ip\{Yt,t) — 1/2. We have checked both numerically and 
analytically that the following picture holds for t and T}, 
large, with a given ratio t/rx larger than 1: 

In the range where x — Yt is of order 1 



■)px{x,t) ~ (j)y{t)ix - Yt), 



(18) 



where v(t) — Yt is the instaneous velocity of the front and 
(j)^ is the solution of the Fisher-KPP equation moving at 
a constant velocity v, i.e. the solution of 



C + V'^'v + (1)1 - (t>v ^ 0, 



with0^(— oo) = and 0„(+cxd) = 1. (The same form ( 18 ) 
is used in 34J.) For definiteness, we normalize such that 
(t)v(0) = 1/2. This determines a unique solution which 
has, if w > 2, the following asymptotics for z — > +oo: 

l-Mz)-B^e-'^"+o{e-'"'), (20) 

where 7 is the smallest solution of 

v = j + -f-\ (21) 

On the other hand, in the range x — Yt ^ 1, 'ip\{x, t) is 
accurately given by the solution of the equation obtained 
by linearizing ([2| around 1: 

l-VA(x,t)~ ^^^^ r e-^"du. (22) 



Then, using the asymptotics of the error function 
/y eicp{—v?)du ~ exp(— X^)/(2X) and requiring that 



( 18|20 l and (22l match in the range 1 < a; - Yt < V^, 

one gets that ^{t) and Yt should satisfy 

i?,(t)e-^W(-^') (23) 

To match the dependence in x—Yt, we need Yt ~ 2t^(t) 
to first order. Then matching the prefactors leads to 



Yt ~ 2^t{t - Tx) - 



Int 

with j{t) 



ln[2V^7(t)i?^(t)] 
lit) 



t-Tx 



t 



(24) 
(25) 



Note that the relation (21 1 is indeed satisfied to leading 
order as Yt ~ j{t) + ^(t)^^. In figure [2] we see that 
the agreement between the leading term in ( |24[ ) and the 
position obtained by integrating numerically with the 
initial condition ^ is quite good. One could also see in 
( 24 1 a 1/t convergence of the A-dependent delay which is 
consistent with the numerical results of figure [l] 




(19) FIG. 2: The prediction (24 1 to the leading order for the posi 



tion of the front is compared to the position measured by 
integrating ([2| with the initial condition ([sjl. As in (16 1, 
TA = -ln(l- A). 



2. For V = 2, 
Aze^^ for large 



For t/rx large, j{t) — * 1 and v{t) 
the solution of (19) satisfies 1 — (p2{z) 
z dH For V slightly larger than 2, the next term 
in the large z expansion (20 1 is 1 — <j)v{z) — Bj e^^^ + 
C-^e"^!^ -I- o(e~^l'^). For consistency in the limit v 2, 
one has 



A 



2(1-7) 



as 7 



(26) 



7(t) 



SO that B 
comes Yt ~ 2t 
gives ( |13|14|16 1 



AtjTx from ( |25|26| . Thus, ([24]) be- 
- TX - (3/2)lnV+IhrA + 0(1), which 
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One can repeat everything if, instead of starting with 
a single particle at the origin, one starts with K of parti- 
cles at positions yi,. . . ,yK- One simply needs to replace 
ip\{x,t) defined in Q by ni<i<xV'A(a^ - yi,t), with a 
similar change for CJo{x,t) = tpo{x,t). As a result, in 
the long time limit, the delay function /(A), and there- 
fore the distances between the rightmost particles remain 
unchanged. This property is remarkable: whatever the 
positions of the initial particles are (as long as there are 
a finite number of them) the limiting average distances 
and probably the whole limiting measure seen from the 
rightmost particle are the same. 

We can also extend all our calculations to more gen- 
eral branching random walks. For example one may con- 
sider a discrete time case where at each time step, every 
particle splits into K new particles, and the position of 
each new particle is shifted from its parent by a random 
amount e drawn from a given distribution /5(e). Apart 
from a few changes, such as (21 ) which is replaced by 



In 



K / p{e) e^' de 



(27) 



which becomes the solution of 



_ In(l-A) 



, dv 

d'y 



(28) 



everything remains unchanged. In particular ( 16|17 l are 
simply divided by the value 70 of 7 which minimizes the 
expression (27 1 of w. Thus, the asymptotics of both the 



delay (16 1 and the distances (17 1 look universal, up to a 



scale factor 79. 

In the present letter we have seen that the distances be- 
tween the rightmost particles at the frontier of a branch- 
ing random walk have statistical properties ( 10|17 1 which 
can be understood as the delay ( 14|15|16" l of a traveling 
wave. Other properties, such as the correlations of these 
distances or even their whole probability distribution can 
also be understood in terms of the delay of a traveling 
wave. For example if Rn^m{x, y, t) is the probability that 
there are n particles at the right of x and m particles at 
the right of y, on can show that 



{t)d it)) 



dxx I dyy 



d'^ Rn,m{x,y,t) 
dxdy 



while the generating function ^„ A"^"'i?„,m(a::, x+c, t) 
defined as in ^ evolves according to the Fisher-KPP 
equation (|2| with a new initial condition. 

A surprising aspect of the present work is that the 
statistics of the leading particles, in the long time limit, 
do not depend on the positions or on the number of par- 
ticles we start with, as long as there is a finite number of 
them. This means that the limiting measure has the fol- 
lowing stability property: if one takes two realizations of 
the leading particles according to this measure and shifts 



one of them by an arbitrary amount, then the superim- 
position of these two realizations gives a new realization 
of the same measure, up to a translation. 
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